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ABSTRACT 

We discuss atmosphere models of HD209458b in light of the recent day-side flux measurement of 
HD209458b's secondary eclipse by Spitzer-MIPS at 24 /im. In addition, we present a revised secondary 
eclipse IRTF upper limit at 2.2 /im which places a stringent constraint on the adjacent H2O absorption 
band depths. These two measurements are complementary because they are both shaped by H2O ab- 
sorption and because the former is on the Wien tail of the planet's thermal emission spectrum and the 
latter is near the thermal emission peak. A wide range of models fit the observational data, confirming 
our basic understanding of hot Jupiter atmospheric physics. Although a range of models are viable, 
some models at the hot and cold end of the plausible temperature range can be ruled out. One class 
of previously unconsidered hot Jupiter atmospheric models that fit the data are those with C/O >1 (as 
Jupiter may have), which have a significant paucity of H2O compared to solar abundance models with 
C/O = 0.5. The models indicate that HD209458b is in a situation intermediate between pure in situ 
reradiation and very efficient redistribution of heat; one which will require a careful treatment of atmo- 
spheric circulation. We discuss how future wavelength-dependent and phase-dependent observations will 
further constrain the atmospheric circulation regime. In the shorter term, additional planned measure- 
ments for HD209458b, especially Spitzer IRAC photometry, should lift many of the model degeneracies. 
Multiwavelength IR observations constrain the atmospheric structure and circulation properties of hot 
Jupiters and thus open a new chapter in quantitative extrasolar planetology. 
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1. INTRODUCTION 

The star HD209458 is the brightest known {V = 7.5) 
to host a short-period giant planet which transits the face 
of the star (Charbonneau et al. 2000; Henry et al. 2000). 
Recently Deming et al. (2005b) detected the secondary 
eclipse of the planet by the star in this system, at a wave- 
length of 24/im using the MIPS instrument on the Spitzer 
Space Telescope. The measured flux decrement during sec- 
ondary eclipse provides a direct measurement of the plan- 
etary day side thermal emission. Together with the recent 
Spitzer IRAC secondary eclipse detection of TrES-1 (Char- 
bonneau et al. 2005), this marks the first measurement of 
the thermal emission from a known extrasolar planet. 

In this paper we examine what constraints can be placed 
on the planetary properties by this historic measurement. 
While our principal constraint is the measured 24 mi- 
cron flux [Fiifaxi) there are several additional observational 
data points useful for the study of the HD209458b atmo- 
sphere. The most relevant additional constraint for the 
planetary day side flux is an upper limit on the H 2 ab- 
sorption band depths on either side of the 2.2 /im contin- 
uum flux peak (AF^^m), reported by Richardson et al. 
(2003). Two data points from transmission spectra of the 



lower atmosphere — which probe the planetary limb — are 
relevant: the Na resonance doublet detection at 0.584 /im 
(Charbonneau et al. 2002) and an upper limit on CO at 
2 /im (Deming et al. 2005a). With these data, we can 
begin to constrain atmosphere models of the extrasolar 
planet HD209458b. 

The new Spitzer data is consistent with a circular orbit 
for HD209458b, based on the equal time spacing between 
primary and secondary eclipses (Deming et al. 2005b). In 
addition, the latest radial velocity determination of the 
eccentricity is essentially zero: e = 0.013 ± 0.009 (Laugh- 
lin et al. 2005). With a circular orbit, the anomalously 
large radius of HD209458b is not explained by an in- 
teracting companion (Bodenheimer et al. 2003; Deming 
et al. 2005b). While atmospheric processes have also been 
hypothesized to play a role in the radius evolution of 
HD209458b (Guillot & Showman 2002), to first-order, the 
atmospheric circulation regime on the ~ 3-day period hot 
Jupiters with radii closer to Jupiter's is not expected to 
differ much (Menou et al. 2003). The -F^m also indicates 
that the planet HD209458b is in synchronous rotation with 
its orbit, presenting the same face to the star at all times. 
Because the timescale for synchronization is much smaller 
than the timescale for orbital circularization (Goldreich & 
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and 9.78 x 1CT 6 ± 1.52 x 1(T 4 respectively. We aver- 
age these two values to obtain an average band depth of 
1.5 x 10~ 5 ±9.8 x 10~ 5 . Since the result is essentially zero, 
we take the error on the mean ~ 1 x 10~ 4 to be the upper 
limit on the band depth on either side of the 2.2 /im flux 
peak at the la level. 

We convert this contrast ratio limit into a flux limit of 
Af2.2/jm = 200yuJy (see Figure 1), using a Kurucz model 
atmosphere for the star (Kurucz 1992) and the HD209458b 
planet-star area ratio (using R p — 1.347 ± 0.060i?j and 
R* = 1.146±O.O5Oi? from Brown et al. 2001). To summa- 
rize, the AF 2 . 2^111 constraint means that the band depths 
from the continuum peak to the absorption trough must 
be < 200/iJy. This proves to be a useful constraint and is 
discussed further in relation to model atmospheres in §3. 

3. hd209458b atmosphere models 

Model atmosphere computations are required for an in- 
terpretation of spectral data because the planetary fluxes 
at different wavelengths are shaped by a variety of physical 
processes. Measurements of the flux level at only one or 
two wavelengths can be very misleading if not interpreted 
in the appropriate model context. 

3.1. Background and Caveats 

Model spectra were computed with a ID plane parallel 
radiative transfer code. This code solves three equations: 
the equation of radiative transfer, the equation of radiative 
equilibrium, and hydrostatic equilibrium to derive three 
unknowns: temperature (as a function of altitude), pres- 
sure (as a function of altitude) and the radiation field (as 
a function of altitude and wavelength). The boundary 
conditions are the stellar radiation at the top of the atmo- 
sphere, and the interior entropy at the bottom of the atmo- 
sphere. The model outputs are temperature, pressure, and 
radiation field (from which the emergent flux can be com- 
puted) . From the wavelength-integrated thermal emission 
flux (F), T eq can be calculated by (F/api) 1 ^ 4 where or 
is the radiation constant, and the geometric albedo 1 can 
be computed by ratioing the visible-wavelength emergent 
flux to the incident stellar flux. 

Although the physics of this ID model is straightfor- 
ward, it is necessary to choose several input parameters, 
which thus constitute uncertainties with various levels of 
impact. The interior entropy is unknown, although it can 
be estimated from evolutionary models. Opacities govern 
the absorption and reemission of radiation and so drive 
the entire radiative transfer problem, determining the ver- 
tical temperature-pressure profile and the planetary flux. 
The choices of metallicity, which atomic and molecular 
species to include, whether to include only equilibrium or 
non-equilibrium chemistry as well, all affect the opacities 
and are thus input parameters. Clouds are the most se- 
rious uncertainty in terms of opacities. They form at a 
the saturation vapor pressure and (for solar system plan- 
ets) can extend 1 to 2 pressure scale heights above this 
"cloud base". Dozens of high temperature condensates 
could potentially exist; only a few are suspected to have a 

1 The Bond albedo (the total radiation reflected in all directions compared to the total incident radiation from the star) cannot be computed 
in these ID models, and is difficult to derive from the computed geometric albedo (reflected radiation relative to that from a flat Lambertian 
surface of the same cross-sectional area as the planet). On most solar system planets the Bond albedo ranges from 10% to 40% less than the 
geometric albedo. See Chamberlain & Huntcn (1987) for more quantitative albedo definitions. 



Soter 1966), the tidally-locked assumption is likely valid. 
The notion of a permanent day side has important con- 
sequences for model atmospheres and data interpretation 
because it means only one side of the planet is being heated 
by the star. 

We begin by describing the HD209458b measure- 
ment and AF2.2fim upper limit in §2. In §3 we describe the 
radiative transfer models and their uncertainties. In §4 we 
present a comparison of the models with the HD209458b 
data as well as a discussion of the relation between the 
observations of TrES-1 and HD209458b. In §5 we con- 
sider the broader implications and near-term prospects for 
both theory (§5.1) and observation (§5.2). We conclude 
in §6 by describing how specific upcoming observations 
of HD209458b will enable a more definitive atmospheric 
characterization. 

2. HD209458B secondary eclipse data 

In this paper we focus on the thermal emission data. 
Both i<24 Mm and AF 2 .2^m are secondary eclipse measure- 
ments which probe the planetary day side thermal emis- 
sion. 

2.1. The Spitzer 2^jim Flux 

The HD209458b secondary eclipse measured by Spitzer 
yields a planet/star flux ratio F 2 4 Mm /F* = 0.00260 ± 
0.00046, a 5.6cr result. During secondary eclipse the planet 
is obscured and only starlight is present. The resulting 
flux decrement yields the ratio of planet to stellar flux. 
Multiplying this by the measured stellar flux at this wave- 
length (21.17 ± 0.11 mJy), the planetary flux is thus mea- 
sured to be i^24pm = 55 ± 10 yuJy. 

2.2. The IRTF 2.2 [im Constraint 

Hot Jupiter atmosphere models generically predict a 
peak in the continuum near 2.2 /im due to the combined in- 
fluence of water and CO absorption (with CO a relatively 
weak component at 2.3 /im). Richardson et al. (2003) 
used the method of "occultation spectroscopy" to search 
for this signature during the HD209458b secondary eclipse 
using the NASA Infrared Telescope Facility (IRTF) SpeX 
(Rayner et al. 2003) spectral data from 1.9 to 4.2 microns 
at a spectral resolution of A A/ A = 1500. We have revisited 
the data to improve the upper limit. 

The data do not provide an absolute upper limit to the 
2.2/xm flux itself. The differential nature of the obser- 
vations means that only a relative measurement among 
the data points in the specified spectral range is possi- 
ble. We therefore calculate the upper limit to the flux 
difference, or band depth, at the location in the spectrum 
where an absorption band is expected. Note that the data 
points are expressed in terms of the star-planet flux ra- 
tio. We use three data bins: a central bin where the 
flux peak is predicted to reside (2.09-2.31/mi), and one 
bin on either side of the flux peak bin where an absorp- 
tion band trough should appear (1. 955-2. 09/im and 2.31- 
2.52/im). The difference between the flux peak bin and 
the absorption trough bins are 2.02 x 10~ 5 ± 1.25 x 10~ 4 
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significant magnitude of opacity. The condensate opacity 
is controlled by particle size distribution and index of re- 
fraction, as well as particle shape. Although cloud models 
have been used by some researchers to compute the parti- 
cle size distribution and cloud vertical extent for spherical 
particles (Ackerman & Marley 2001; Cooper et al. 2003), 
they still have as a free parameter the type of condensate 
and the amount of condensed material in the cloud (e.g., 
by specifying the efficiency of sedimentation). All extra- 
solar planet atmospheres models in the literature consider 
horizontally homogeneous clouds (i.e., covering the entire 
planet). 

Finally, for ID models we must also make an assumption 
of how the incident stellar radiation is redistributed across 
the surface of the planet. We use the parameter / (de- 
scribed in §4.1) as the proxy for atmospheric circulation: 
/ = 1 if the absorbed stellar radiation is redistributed 
evenly throughout the planet's atmosphere (e.g., due to 
strong winds rapidly redistributing the heat) and / = 2 
if only the heated day side reradiates the energy (for the 
latter T eq is for the day side only). The implementation in 
the models is achieved by multiplying the incoming stellar 
radiation 2 by a factor of //4. 

3.2. Model Choices 

In terms of choices for the models presented here, we 
use: a Kurucz model atmosphere for the stellar radia- 
tion (Kurucz 1992); solar abundances (HD209458A is close 
to solar metallicity with [Fe/H] = 0.04 (Gonzalez et al. 
2001); and line opacities of H 2 0, CH4, CO, Na, K, and 
collision-induced absorption opacities of H 2 -H 2 (Borysow 
et al. 2001) and H 2 -He (J0rgensen et al. 2000). We adopt 
a planetary interior energy flux 3 corresponding to an ef- 
fective temperature T e g = 500 K. A Gibbs Free energy 
minimization code was used to compute chemical equilib- 
rium (Seager 1999). H 2 Rayleigh scattering and scattering 
and absorption by MgSiOa and Fe vertically and horizon- 
tally uniform and homogenous clouds are included with 
a specified particle size distribution (a log normal distri- 
bution with particle radius of 2 ^im and a = 0.05) and 
vertical extent (two pressure scale heights above the cloud 
base) are also included. Note that for computational ef- 
ficiency we use a separate, higher spectral resolution, ra- 
diative transfer code with the same input assumptions to 
compute the final emergent spectra based on the T/P pro- 
files. A more detailed description of the code and opacity 
sources can be found in Seager (1999) and Seager et al. 
(2000). (See Barman et al. 2001; Sudarsky et al. 2003; 
Burrows et al. 2005; Fortney et al. 2005, for choices made 
by other modelers). 

With all of the above uncertainties there are many 
choices that lead to a wide range of models for a given 
planet. Rather than an exhaustive search of the model pa- 
rameter space, we consider three models (Figure 1) which 
span a range of reasonable possibilities, to interpret the 
HD209458b F 24Aim and AF 2 . 2Mm data. Model 1 is possibly 
the simplest case one can consider; an / = 2, cloud-free 
model with solar abundances. The second case, model 2, 
is a similarly simple model, but now with / = 1, so that 



some energy is redistributed to the night side. Model 3 
illustrates the effects of thick clouds; it is an / = 1 model 
with thick, high clouds given by the parameters described 
above. 

4. DATA INTERPRETATION 

4.1. Planetary Equilibrium Temperature 

The planetary equilibrium effective temperature T eq is a 
fundamental parameter of the planetary atmosphere and 
is useful as a proxy for the global temperature of a planet 
atmosphere. To first order T eq tells us how much radia- 
tion is being absorbed by the dayside of the planet. In this 
subsection we relate the measured F 2 4 Mm to T eq . 

The measured flux -F 2 4 Mm = 55±10/iJy can be converted 
into a 24/im brightness temperature T 2 4 Mm = 11 30 ±150 K 
(Doming et al. 2005b). The brightness temperature is the 
temperature of a blackbody emitting the flux equivalent 
at a specific wavelength and is not necessarily an accurate 
representation of the planet's equilibrium effective temper- 
ature T eq . Only if the planet radiates as a blackbody is 

The spectrum of a hot Jupiter, however, is expected 
to be very different from a blackbody due primarily to 
strong H 2 absorption throughout near-IR and mid-IR 
wavelengths. In particular, in the mid-IR (which is the 
spectral region probed by MIPS) H 2 has strong continu- 
ous absorption such that there is no true planetary contin- 
uum expected from models at the MIPS wavelength bands, 
even though the relatively flat spectrum makes it appear 
so. In other words, the H 2 vapor absorption in the 24 
region substantially depresses the planetary flux. This 
implies that if H 2 vapor absorption is present, the true 
T eq > 1130 K. Figure 1 shows the theoretical spectrum of 
our three HD209458b example models; each of these mod- 
els match F 24Atm even though their T eq (1700 K, 1420 K, 
and 1450 K respectively) are all hotter than T 2 4 tim . We 
also show the blackbody spectrum for these temperatures 
to demonstrate how significantly is suppressed by 

H 2 absorption in these models (a flux depression equiv- 
alent to a decrease of ~ 300-600 K in T 2 4 lim ). 

The true equilibrium temperature is regulated by sev- 
eral competing physical effects — the amount of irradiation 
received from the host star, the fraction that is simply re- 
flected rather than absorbed (governed by the Bond albedo 
Ab), and the amount of energy circulated away from the 
hotter regions by hydrodynamic processes in the atmo- 
sphere. Considering energy balance between the absorbed 
and reemitted flux for a planet with zero eccentricity and 
ignoring any (likely much smaller) contribution from the 
planetary interior, 

V>2 [f(l-A B )f\ (1) 



T — T 



2a J 



where i?» is the stellar radius, T* is the stellar effective 
temperature, a the semi- major axis, and Ab the Bond 
albedo. Here / is the proxy for atmospheric circulation, 
derived by considering whether the incident stellar radi- 
ation is reemitted into 47r steradians (/ — 1) if the ab- 
sorbed stellar radiation is redistributed evenly throughout 

2 The planet intercepts stellar radiation in a cross sectional area of 7rR 2 and reradiates the energy into either 27riJ 2 (/ = 2) or 47rR 2 (/ = 1). 

3 Although interior models predict an interior T e g ~ 100 K a hotter temperature is required to explain HD209458b's large radius (Guillot & 
Showman 2002) and motivates our adopted planetary interior energy flux corresponding to the T e g of 500 K. 



4 



the planet's atmosphere (e.g., due to strong winds rapidly 
redistributing the heat) or reemitted into 27T steradians 
(/ = 2) if only the heated day side reradiates the energy 
(for the latter T eq is for the day side only). Figure 2a 
shows T eq for the HD209458 parameters (T* = 6000 K, 
R„ = 1.18i? , and a = 0.046 AU; Mazeh et al. 2000; Cody 
& Sasselov 2002). 

An equilibrium temperature T eq = 1700 K is at the hot 
end of the range for the HD209458 system parameters, as 



shown in Figure 2. If our model 1 with T eq = 1700 K were 
the unique, correct planetary model, this hot T eq would 
have significant implications: it is only plausible with both 
a low Ab and no heat redistribution (i.e., / = 2). At the 
low end of the plausible temperature range, T eq = 1130 K 
is unlikely. A T eq = 1130 K model would require that 
HD209458b has a very high Bond albedo (> 0.65), and 
would probably require / = 1 because Ab > 0.8 is dif- 
ficult to attain. Most of the solar system planets have 
Ab < 0.35, with the exceptions of Venus (Ab = 0.75) and 
Pluto (A B = 0.4) (de Pater & Lissauer 2001). 

4.2. Constraints on Models 

Model 1 (solar metallicity, no clouds, no energy redistri- 
bution) is characterized by strong near-IR H2O vapor ab- 
sorption bands. The H 2 absorption also causes a strong 
depression of i^/im- With the strong alkali metal absorp- 
tion dominating H2 Rayleigh scattering at visible wave- 
lengths, this model has a low albedo (geometric albedo = 
0.15) and hence high T eq = 1700 K. 

The Na detection and CO upper limit transmission spec- 
troscopy are both too weak to match a solar abundance 
cloud-free chemical equilibrium model (Charbonneau et al. 
2002; Doming et al. 2005b) and therefore in general pro- 
vide a useful model constraint. In this cloudless case the 
Na could be transported to the night side and converted 
to Na2S at the lower temperatures (Guillot & Showman 
2002; Iro et al. 2005); similarly CO would be converted to 
CH4 on the night side, thus also potentially satisfying the 
limits from transmission spectroscopy. It is important to 
keep in mind that the transmission spectroscopy probes 
the planetary limb; with / = 2, the day and night side are 
assumed to be different temperatures and so the planetary 
limb could have different conditions from those computed 
on the day side, for example a high cloud at the limb due 
to the falling temperatures away from the substellar point. 
This high cloud would also satisfy the Na and CO limits 
from transmission spectroscopy. 

We use this T eq = 1700 K model as an example of a 
model that is ruled out by the data. Although this model 
satisfies i^^m an d can satisfy the transmission spectra 
data, it is ruled out (at the 5cr level) by the AF2.2^m 
data. The strong water vapor absorption bands adja- 
cent to 2.2/im are in conflict with the Ai^^m limit. 
The A_F2.2/jm constraint means that the band depths from 
the continnuum peak to the absorption trough must be 
< 200/iJy, whereas in this model the band depths are 
~ lOOO^Jy. We used this simple model in Deming et al. 
(2005b) primarily to illustrate that T eq does not have to 

4 Recall that in LTE with the source function equal to a Planck function, the solution to the ID plane parallel radiative transfer equation is 
F(u) = 2tt Jq 1 7(r, /i, v)fid/i = 2ir J °° B(T(t), v) cxp(— r/fi)drdfi, where T is a function of the altitude expressed on an optical depth scale r, 
fi = cos 9 where 8 is the angle away from the normal, I is the intensity, and v is frequency. Taking B(T(t), v) as a constant out of the integral 
we see that for an isothermal atmosphere F(y) = ttB(T, v). 



equal T24 flm . The fact that model 1 is ruled out by our 
revised 2.2 /<m upper limit shows that AF 2 .2 Mm is a useful 
constraint on models. 

Model 2 (in which energy is now redistributed to the 
night side) has a lower computed T eq = 1420 K. While this 
model fits F 2 4 Mm , it only fits the Ai*2.2/tm constraint at the 
3cr level due to the strong H 2 absorption features. In the 
context of this / = 1 cloud-free, solar abundance model 
we note that the atmospheric circulation could allow for 
substantially varying temperatures at different horizontal 
locations (e.g., Showman & Guillot 2002; Cho et al. 2003; 
Cooper & Showman 2005) and so more detailed models 
may show partial cloud cover and lower Na abundances 
from molecular condensation (which may also better fit 
the weak Na and CO transmission spectra) . 

In Model 3 we allow for formation of thick, absorp- 
tive clouds and so, in contrast to models 1 and 2, it 
has very weak absorption band features with a computed 
T eq = 1450 K. The absorption and remission and scatter- 
ing from largely grey cloud particles causes the absorption 
bands, including the continuous H 2 absorption at 24 mi- 
crons, to be much shallower than in a cloud-free model. 
The cloud bases are computed to be at 5 and 10 millibar 
for MgSiOs and Fe respectively, with the cloud tops two 
pressure scale heights above at 0.8 and 1.3 millibar; this is 
consistent with both the lower-than-expected Na and CO 
transmission spectra. As seen from Figure 1 this model fits 
both the Ai^jum and F24 fl m data. We emphasize that the 
amount of condensates and their opacity strength (con- 
trolled by particle type and size distribution for spherical 
particles) hugely affects the overall spectrum; the conden- 
sate opacity is competing with the H2O opacity to de- 
termine the water absorption band strengths. Different 
choices for the particle type, fraction of species condensed, 
cloud vertical scale height, and particle size distribution 
will all affect the overall opacity. Just one example is that 
larger particle sizes are likely to be more reflective and 
would not weaken the absorption bands nearly as much as 
absorptive particles. 

4.3. Consequences of Model Constraints 

Based on our investigation of models 1 through 3, we can 
draw some general conclusions. We showed above that the 
T eq — 1700 K model, under our adopted input assumptions 
(described in §4), is ruled out based on the AF; 



2.2/im 



ab- 



sorption band depth constraint. Other hot models are also 
ruled out by the data. First, a truly isothermal cloud- free 
atmosphere would satisfy the Ai^^m because a vertically 
isothermal LTE atmosphere does not have any spectral 
features — i.e., it is a blackbody spectrum 4 . However, such 
an isothermal T eq = 1700 K model is ruled out by the 
-F24/xm data even considering the 3a error on T 2 4 Aim which 
permits a value as high as 1580 K. Second, most hot mod- 
els in between the isothermal atmosphere and the strong 
vertical temperature gradient atmosphere are likely ruled 
out; as the vertical temperature gradient gets shallower 
and satisfies AF 2 .2^m with a weaker absorption band, the 
24 /im model flux gets higher and away from the mea- 
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sured i^^m- Third, a hot atmosphere with clouds that are 
highly absorptive (low Ab), such as Fe clouds, could also 
result in a blackbody spectrum (even without an isother- 
mal atmosphere) — again a hot blackbody or near black- 
body spectrum is ruled out because T 2 4^ m = 1130 K is 
lower than resulting values of T eq . 

A second category of models that are difficult to fit are 
those with T eq at the cold end of the plausible range (equa- 
tion (1) and Figure 2). As described above, for HD209458b 
to have a T eq near 1130K, a very high Ab > 0.65 is re- 
quired. Such a high As could be possible only from a 
narrow range of cloud parameters and is not attainable 
in the absence of clouds — i.e., due to Rayleigh scattering 
alone. 

Both of our / = 1 and / = 2 cloud-free models have a 
strong H 2 absorption band on either side of the 2.2 ztm 
flux peak that do not fit the Ai*2.2^m constraint at the la 
level. What could cause weaker H2O absorption features 
to fit Ai^^/^m? Absorption and remission by cloud parti- 
cles weaken the 2.2 /im absorption feature and are hence 
consistent with the data. A more isothermal atmosphere 
would also have weaker absorption features (see above). 
Decreasing the water vapor abundance could also serve 
to weaken the water vapor absorption features adjacent to 
2.2 ztm. This situation is hard to physically realize in a so- 
lar abundance atmosphere. Photodissociation from stellar 
UV radiation affects water only in the upper atmosphere, 
close to the microbar pressure level (Liang et al. 2003) 
where the low H 2 abundance would have little effect on 
the overall planetary spectrum. 

4.4. Low H2O Abundance in the C/O > 1 Regime 

A possible scenario with a very low abundance of water 
vapor, that would satisfy Ai*2.2f*m) is one with atmospheric 
abundances of the carbon-to-oxygen ratio C/O > 1. This 
is in contrast to the solar C/O ratio used in all extrasolar 
planet models published so far (0.5 or 0.42, Allende Pri- 
eto et al. 2002; Anders & Grevesse 1989, respectively). In 
addition to being metal-rich overall, Jupiter is likely en- 
riched in C with C/O = 1.8 (Loddcrs 2004). In addition, 
Saturn has less oxygen enrichment than other heavy el- 
ements, including carbon (Visscher & Fegley 2004). The 
enhanced metallicity on Jupiter and the other solar system 
planets is thought to come from post-formation planetesi- 
mal accretion and may well be applicable to hot Jupiters; 
enriched C compared to O over solar abundances would 
require pollution from carbon-rich planetesimals. 

A high C/O ratio > 1 can have a large effect on the 
spectral signatures of hot planetary atmospheres (see Fig- 
ure 3 and Figure 1 in Kuchner & Seager 2005, Seager & 
Kuchner, in preparation). This is because at high temper- 
atures the dominant carbon-bearing molecule is CO, and 
in equilibrium formation of CO is chemically favored in- 
stead of H2O. In an extreme case, then, with all of the 
O tied up in CO, there will be no H2O vapor present. 
In reality, depending on the temperature-pressure profile, 
there will be a small amount of water vapor present, but a 
few to > eight orders of magnitude below the CO number 
density. This extreme paucity of H 2 is significant, be- 

We chose the value C/O = 1.01 because graphite clouds tend to 
1997). 

6 Recall that T eq is for the planetary day side. 



cause the hot Jupiter spectra are normally predominantly 
shaped by H 2 absorption bands. In the C/O > 1 case 
the planetary spectrum will instead be dominated by CO 
absorption and by CH4 absorption (which forms at higher 
temperatures than for C/O < 1). (Note that for 0.5 < 
C/O <1, the water vapor will decrease much less signifi- 
cantly, by a factor of a few to 100 respectively. See Fortney 
et al. (2005) for a discussion of the effects of C/O = 0.7.) 
In contrast to the hot Jupiter case, a C/O ratio > 1 has 
a much less significant effect on the spectral signatures of 
planetary atmospheres cooler than the hot Jupiters (below 
T eq ~ 1000 K). In this case CH4 is the dominant carbon- 
bearing molecule regardless of the C/O value; whatever O 
is present is free to form H 2 0. 

Any hot Jupiters with C/O > 1 should have a very low 
amount of water vapor, more so for those at the hot T eq 
range and those with hot upper atmospheres. Figure 3 
shows an estimated thermal emission spectrum (using a 
scaled model 1 T/P profile to reproduce a T eq = 1600 K) 
with (condensate-free) chemical equilibrium computed for 
a value 5 of C/O = 1.01. More careful modeling is needed to 
explore the consequences of a carbon-enriched hot Jupiter 
atmosphere, particularly for condensate formation which 
is affected by the O abundance. For C/O > 1 graphite 
and SiC clouds can form which are not present in C/O 
< 1 models. 

4.5. HD209458b and TReS-1 Comparison 

We now turn briefly to a discussion of the thermal 
emission detected from TrES-1, and how it relates to 
HD209458b day side thermal emission. The secondary 
eclipse detection of TrES-1 (Alonso et al. 2004) by Char- 
bonneau et al. (2005) at Spitzer IRAC bands (4.5 and 8 
/zm) was reported at the same time as the HD209458b sec- 
ondary eclipse detection. A planetary brightness tempera- 
ture was reported at each bandpass: T 4 . 5A , m = 1010 ±60 K 
and Tg[i m = 1230 ± 110 K, with an average Tttes-i = 
1060 ± 60 K. In this section we consider the TrES-1 data 
as it relates to HD209458b. 

As for HD209458b, Tti-es-i does not constrain Ab due 
to the unknown atmospheric redistribution of absorbed 
stellar radiation (/ in equation (1); c.f. Charbonneau et al. 
2005). If T TrE s-i = Teq = 1060 K (shown in Figure 2b 
for the TrES-1 parameters R* = 0.85i? Q , T* = 5250 K, 
a = 0.039 AU; Alonso et al. 2004), then 0.3 < A B < 0.7, a 
large range. Using the same arguments as for HD209458b 
(non-blackbody-nature of the spectrum due primarily to 
H 2 absorption) the T eq is likely to be hotter than HOOK. 
A very interesting consequence of TrES-1 stellar temper- 
ature and planet semi-major axis is that if T eq > 12QQK 
an f = 1 model is ruled out e , otherwise we would have the 
unphysical situation Ab < 0. 

We emphasize that HD209458b and TrES-1 arc likely 
to be different planets in terms of their atmospheric struc- 
ture. We caution that the similar measured brightness 
temperatures of TrES-1 and HD209458b are misleading; 
they are at different wavelengths and probably probe dif- 
ferent atmospheric altitudes/temperatures. Furthermore, 
TrES-1 and HD209458b are theoretically likely to have 

te the gaseous carbon close to a C/O ratio of 1 (Loddcrs & Fegley 
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different T eq due to the former's intrinsically fainter par- 
ent star (KOV (T cff = 5250 K) vs. GOV (T cS = 6000 K) 
respectively; see Figure 2). This temperature difference 
could have important consequences in terms of cloud alti- 
tude and the CO/CH4 abundance ratio. A cold model for 
TrES-1 would have more CH 4 in the atmosphere than the 
hotter HD209458b models, for the same abundances. The 
two planets are also different in their radii (for HD209458b 
R P = 1.347±0.06ORj (Brown et al. 2001) and for TrES-1 
R P = IMtl^Rj (Alonso et al. 2004)). 

Although the two planets are different, it is still instruc- 
tive to compare the TrES-1 data to the HD209458b mod- 
els (Figure 1), since the solar abundance models have the 
same generic H2O and CO absorption features at 4.5 and 
8 /xm. It is certainly puzzling that the 8/jm flux is > 
the 4.5^tm flux, because one generally expects more ther- 
mal emission at shorter wavelengths (see Figure 1). The 
TrES-1 data constrain the relative flux ratios of a shallow 
H 2 band at 8ttm and a deeper H 2 + CO band at 4.5 
microns. For most model fluxes that match the 8^m band 
the model fluxes are too high in the 4.5/im band (Char- 
bonneau et al. 2005) but only marginally so considering 
the 3<7 errors on the data. 

One solution to the IRAC 8/im and 4.5/im data dis- 
crepancy (Charbonneau et al. 2005) is a missing source 
of opacity at 4.4/ttm to cause a deeper absorption band. 
This discrepancy is at face value opposite to the require- 
ment of weak H 2 bands for HD209458b from AF 2 . 2Aim . 
Additional H 2 opacity would not likely solve the IRAC 
measurement discrepancy, because H 2 is present in both 
observed IRAC bands. An increased CO abundance with- 
out an increased H 2 abundance could potentially solve 
the problem because CO absorbs in the 4.4 /jm band and 
not the 8 /xm. A higher CO number density and lower 
H 2 number density could be realized with a higher C/O 
ratio than solar (but not necessarily C/O >1), as long as 
not too much CH4 is present (which also absorbs in the 
8/zm IRAC band and whose number density grows as C/O 
increases) (Seager & Kuchner, in preparation). See Fort- 
ney et al. (2005) for a discussion of metallicity and varying 
C/O ratios for TReS-1 models. A second solution to fit the 
nearly equal fluxes in the IRAC 8/im and 4.5/xm bands is 
a close-to- vertically isothermal model with weaker absorp- 
tion bands. Such a model would have to be a cool model to 
fit the data, as illustrated by the blackbody fit in Figure 1 
and in Charbonneau et al. (2005). We note that a cooler 
model fit indicates that efficient redistribution of absorbed 
stellar radiation (i.e., / — 1) is likely to be present whereas 
for a hotter model it is not (see above) . It is therefore im- 
portant for future observations to distinguish between the 
best fit models: a hotter model with strong absorption 
bands or a colder model with weak absorption bands. 

5. ATMOSPHERIC CIRCULATION: SELF-CONSISTENT 
MODELS AND FUTURE OBSERVATIONAL DIAGNOSTICS 

We now turn from interpreting data, to examining fu- 
ture prospects for theoretical models. One of the major 
simplifications in the above models is the basic assumption 
about how the absorbed stellar radiation is redistributed 

7 On the planetary dayside, vertical heat transport by convection is not 
imposes a shallow vertical temperature gradient. Additionally, we note 
near the photosphere (i.e., t ~ 1; e.g., Iro et al. 2005). 



throughout the planet atmosphere, with an extreme / = 1 
or / = 2 assumption. Indeed, this is one of the most com- 
pelling issues to understanding the hot Jupiter data and 
models alike. With a permanent day side facing the star 
from tidal locking and intense stellar radiation from the 
very small semi-major axes (< 0.05 AU) strong winds are 
expected to develop to advect the heat away from the plan- 
etary day side (e.g., Showman & Guillot 2002; Cho et al. 
2003; Menou et al. 2003; Cooper & Showman 2005). In this 
section we describe why such self-consistent coupled radia- 
tive transfer/atmospheric circulation models are required, 
and discuss the possibility of constraining the atmospheric 
circulation regime on HD209458b (and other hot Jupiters 
by extension) from wavelength- and time-dependent pho- 
tosphcric signatures such as F 2 4 fim . 

5.1. Radiative and Advective Regimes 

In the simplest terms, the influence of atmospheric 
circulation on the planetary atmosphere horizontal 
temperature-pressure structure can be discussed as a com- 
petition between a typical radiative timescale (r ra d) and a 
typical advective timescale (r a( j v ) (e.g., Showman & Guil- 
lot 2002; Iro et al. 2005). One can roughly estimate the 
relevant timescales from the following relations 7 : 

Trad ~7i^ (2) 

Tadv ~ — , (3) 

where P is the pressure, g is the surface gravity, cp is the 
heat capacity, and U is the characteristic wind speed (a 
priori unknown). We are interested in the ratio of these 
quantities at the location where most of the incoming en- 
ergy is deposited. This is determined largely by the opac- 
ity of the atmosphere at visible wavelengths. 

One can ask: is the bulk of the stellar radiation absorbed 
in the radiative regime r ra d/T a d v -C 1, or in the advective 
regime T ra d/r a d v 3> 1? If the atmosphere is in the radiative 
regime, the bulk of the absorbed energy would be reemit- 
ted before being advected to the night side. In this case 
(which corresponds to our / = 2 model) we would expect 
a strongly phase-dependent horizontal temperature gradi- 
ent centered on the substellar point, and a strong day- 
night temperature difference. If, on the contrary, the bulk 
of the stellar energy is absorbed in the advective regime, 
then atmospheric circulation (Showman & Guillot 2002; 
Cho et al. 2003; Cooper & Showman 2005) is directly re- 
quired to understand the horizontal heat transfer. The 
horizontal temperature field should then be much more 
uniform (corresponding to our / = 1 model) because heat 
is transported and redistributed efficiently over the entire 
planet. 

The main issue in understanding in which regime we 
are is the unknown wind speed U. If we consider U = 
1000 m/s throughout the atmosphere (for illustration pur- 
poses) we can investigate the regime where the bulk of 
stellar radiation is absorbed in our models. For our / = 2 
cloud- free model 1, (Figure 4a) most of the stellar radia- 
tion is absorbed at or above 1 bar. This is right on the 

important until very deep in the atmosphere because stellar irradiation 
that the definition of the radiative timescale adopted here is valid only 
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boundary near T rac i/T ac i v ~ 1- For a smaller U , we would 
be in the radiative regime, and for a larger U in the ad- 
vective regime (since r ra d/T a d v oc U). On the other hand, 
in our / = 1 model 3 with very thick clouds, much of 
the stellar radiation is absorbed much higher in the atmo- 
sphere, around 5 mbar (Figure 4a). In this case, the stel- 
lar energy is absorbed well into the radiative regime. One 
notices immediately that this conclusion contradicts the 
basic assumption of the / = 1 model, namely that energy 
is efficiently redistributed! This ambiguity is probably not 
fatal since the model improvements discussed below will 
alleviate this. 

Considering windspeeds different from U = 1000 m/s we 
see that, according to the simple timescale comparison, ef- 
ficient transport of heat in atmospheric regions with very 
small r ra d requires a supersonic speed (Cooper & Showman 
2005) 8 . While there is no a priori reason to disregard su- 
personic winds from first principles, such a situation would 
be remarkable because all observed wind speeds on solar 
system planets are well in the subsonic regime. This un- 
usual situation could be attributed to the extremely close 
proximity to the star and correspondingly higher temper- 
atures and shorter T rac [ values for hot Jupiters (see equa- 
tion (2)) 9 From a theoretical viewpoint, if circulation with 
large-scale supersonic winds exist, it poses a very inter- 
esting problem for flow adjustment dynamics. Supersonic 
winds are not required, however, to redistribute the en- 
ergy in the short T ra d regime; other mechanisms such as 
subsonic winds/advection, atmospheric waves, or horizon- 
tal radiative transfer can cause energy to be transported 
or to diffuse away from the hot planetary day side. 

All of the above indicate the need for coupled radiative- 
atmospheric circulation models. In particular we need 
to address the unknown value of U, the fact that the 
bulk of the stellar radiation is absorbed near the bound- 
ary between the radiative and advective regimes and 
the apparent contradictions in ID model assumptions 
for heat redistribution. Furthermore, the above simple 
TVad/Tadv timescale argument cannot replace combined ra- 
diative transfer-atmospheric circulation calculations be- 
cause these two atmospheric timescales are actually non- 
linear ly coupled: atmospheric winds are driven by pressure 
gradients due to uneven radiative heating in this context of 
strong insolation. The value of U in equation (3) depends 
on T (wind forcing) while T itself depends on U because 
of heat advection by the winds. 

What do we mean by a coupled radiative trans- 
fer/atmospheric circulation model? The ID emergent 
spectra (§4) depend on the vertical temperature gradi- 
ent. This temperature gradient, in fact, depends on the 
atmospheric circulation to determine the temperatures at 
various altitudes. Further, in the atmospheric circulation 
picture a ID model is not valid, and a computation of 
different vertical (i.e., radial) T/P profiles away from the 
substellar point, followed by the subsequent hemispherical 
integration of the emergent spectra is necessary. As such, 
no truly self-consistent models yet exist in the literature. 
As more data accrues, such consistent models will become 
necessary. 



5.2. Wavelength-Dependent, Phase-Dependent 
Observational Diagnostics 

Infrared observations of hot Jupiters can provide data 
to help constrain the atmospheric circulation via mea- 
sured horizontal temperature gradients. These observa- 
tions include: day-night temperature differences measured 
for transiting hot Jupiters; IR phase curves as a hot Jupiter 
planet orbits its parent star; and the rate of secondary 
eclipse ingress and egress which will differ if the planetary 
dayside is nonuniform in temperature. 

Assuming that we can distinguish observationally be- 
tween strong horizontal temperature gradients and a uni- 
form horizontal temperature, what can be told about at- 
mospheric circulation? In this section we describe how 
wavelength-dependent, phase-dependent observations can 
potentially be used to constrain the planet's atmospheric 
circulation regime. 

The key point is that the altitude of thermal emission 
is wavelength dependent, due to wavelength-dependent 
opacities. This means that at different wavelengths we 
can potentially "see" altitudes with different ratios of the 
radiative to advective timescales. Additionally, the stel- 
lar radiation is absorbed at visible wavelengths, heats 
the atmosphere, and is reemitted at infrared wavelengths. 
The altitude where visible radiation is absorbed compared 
to the altitude where thermal reemission occurs are only 
weakly coupled in the absence of clouds: the chemistry 
of molecular and atomic species (Na, K, H 2 -H 2 collision- 
induced opacities) which control the absorption of stellar 
irradiation at visible wavelengths is not directly coupled 
to the chemistry of molecules that control thermal ree- 
mission at IR wavelengths (predominantly H 2 0). The IR 
spectrum, therefore, allows us to probe different atmo- 
spheric heights, independently of the altitude where the 
bulk of stellar radiation is absorbed. The altitude depen- 
dence of thermal emission is shown explicitly in Figure 4a 
for the models previously considered. Thermal emission 
at 24 fim occurs relatively high in the planetary atmo- 
sphere (P ~ 10 mbar). Wavelengths at the continuum be- 
tween H 2 bands (in particular, in the 1-2/xm range and 
at the 4/xm continuum peak) probe much deeper altitudes 
closer to 1 bar. Thus, we could potentially probe different 
regimes of atmospheric circulation in a single atmosphere 
using multi-wavelength observations. 

We consider two extreme cases, where the radiation is 
absorbed deep in the atmosphere (~ 1 bar) or high in the 
atmosphere (near millibar pressures). 

Building on our discussion in the previous section, if the 
stellar radiation is absorbed deep in the planetary atmo- 
sphere, at pressures > 1 bar, we are likely to be in the 
advective regime. In this regime, redistribution of the ab- 
sorbed stellar energy by the atmospheric circulation will 
be significant. Higher up in the atmosphere, the tempera- 
ture field will be determined by a combination of horizontal 
advection, vertical radiative transport through the various 
atmospheric layers, and horizontal and vertical wave prop- 
agation, but the initial redistribution of heat at the absorp- 
tion levels should have a significant impact on the overall 
temperature structure in the entire atmosphere. Different 



8 The sound speed is ~ 2, 500 m/s at ~ 1500 K. 

9 One should exercise caution, however, since Jupiter, the giant planet with the weakest cloud-deck winds in our solar system, is also the closest 
to the Sun (See Table 1 in Menou et al. 2003). 
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predictions for the temperature field expected in this ad- 
vective regime exist (Showman & Guillot 2002; Cooper & 
Showman 2005; Cho et al. 2003, Cho et al., in preparation) 
and they should ultimately be testable with observations 
in the future once they are combined with detailed spectral 
predictions. 

If the stellar radiation is absorbed at high atmospheric 
altitudes, in the radiative regime, different observational 
signatures are expected. In this case the energy is expected 
to be recmittcd before it can be advected to the planetary 
night side. One expects a strong dayside horizontal tem- 
perature gradient, centered on the substellar point. Such 
a strong temperature gradient could be detected from the 
rate of ingress (and egress) which would be different from 
the uniform horizontal temperature case. An additional 
observational diagnostic is a very strong day-night tem- 
perature difference, which may be detectable by accurate 
flux differencing from near primary and secondary eclipses. 

If we believe the stellar irradiation is absorbed very high 
in the atmosphere where the radiative timescale is short, 
and no strong temperature gradients are inferred (particu- 
larly none centered on the substellar point) then an inter- 
esting possibility exists: very high winds which may even 
be supersonic. However, other processes may also mitigate 
the extreme temperature gradient (see §5.1). An observa- 
tional diagnostic to differentiate between supersonic winds 
and subsonic processes is Doppler broadening of atomic 
lines which could potentially be detected in future observa- 
tions of primary transit transmission spectral lines (Brown 
et al. 2001). 

In short, the prospect of multiple wavelength, phase- 
dependent observations offers the possibility of building a 
picture of the three dimensional temperature field on the 
surface of planets such as HD209458b and of constraining 
the nature of the atmospheric circulation. 

6. SUMMARY AND NEAR FUTURE PROSPECTS 

We have shown that, to first order, ID radiative transfer 
models can fit the HD209458b thermal emission data at 24 
microns (see also Burrows et al. 2004, 2005; Fortney et al. 
2005). This confirms not only that the models work, but 
also that hot Jupiters are indeed hot and likely heated ex- 
ternally by their parent stars. If we conservatively adopt 
the 3<T error bars (Figure 1), a wide range of models fit 
the data and we are far from a unique interpretation of 
the atmosphere of HD209458b. We have described a pre- 
viously unconsidered set of models that can fit the data, 
hot Jupiter atmosphere models with C/O >1 which have a 
significantly different chemical equilibrium than the C/O 
= 0.5 solar abundance models, in particular an extremely 
low abundance of water vapor. Despite the viability of a 
wide range of models the data do allow us to begin to rule 



out some models, specifically hot, cloud-free models and 
models at the cold end of the plausible T eq range. Because 
the very hottest models are ruled out — models only pos- 
sible with / = 2 — a situation intermediate between pure 
radiative equilibrium and very efficient redistribution of 
heat is likely, which will require a careful treatment of 
radiative-advective processes. 

Observations planned for the next year will help towards 
a more definitive characterization of HD209458b. These 
observations are listed in Table 1. The most useful planned 
measurements will be the Spitzcr IRAC 3.6 /jm and 4.5 /im 
bands. From Figure 1 we see that the 3.6 /zm band will 
quantify the overall flux level in a continuum window. The 
IRAC 4.5 /im band will quantify the flux difference in the 
CO and H2O absorption feature, and should support the 
AF 2 .2 M m constraint of a weak H 2 band. Because we be- 
lieve that the very hottest models are ruled out, a geo- 
metric albedo constraint from the visible-wavelength sec- 
ondary eclipse will be complementary in ruling out cold 
models, especially a limit < 0.5 (by MOST or HST; see 
Table 1). A little further into the future, SOFIA could po- 
tentially measure secondary eclipse closer to the planet's 
thermal emission peak at < 4 ^m which will help measure 
the T eq . 

There are many compelling questions concerning the hot 
Jupiters. Do their atmospheres have the basic compo- 
sition we are assuming, such as abundant water vapor? 
Have processes such as atmospheric escape of light gases 
or non-equilibrium chemistry affected the atmosphere in 
ways that have not yet been considered? Is the planet 
mctallicity or C/O ratio greater than solar and an indi- 
cator of planet formation conditions? What is the atmo- 
spheric circulation like on these planets which exist in a 
radiation forcing regime unlike any planets in our solar 
system — can stable supersonic winds be present? With 
the upcoming data on HD209458b and coupled radiative- 
transfer/atmospheric models we can begin to answer these 
questions. 
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Instrument 


Wavelength (^m) 


Diagnostic 


Date 


Comments 


Ref 


Spitzcr IRAC 


3.6, 4.5, 8, 10 


flux level, H 2 + CO 


2005 or 06 


secondary eclipse 


1 


Spitzcr MIPS 


24 


night side T or constraint 


2005 


primary eclipse 


2 


MOST 


0.3-0.8 


single band photometry 
geometric albedo or upper limit 


2004 


data in hand 


3 


HST STIS 


UV-near IR 


spectrophotometry 

geometric albedo or upper limit 


2003 


data in hand 


4 


HST STIS 


UV-near IR 


transmission spectra 

Rayleigh scattering, Na, K, H 2 


2001 


data in hand 


5 


HST NICMOS 


1-2 


transmission spectra 
H 2 


2004 


data in hand 


6 



Table 1 

Upcoming observational measurements for the HD209458b atmosphere. References: (1) Spitzer IRAC GTO 
program PI G. Fazio; (2) Spitzer GO program 3405 PI S. Seager; (3) Walker et al. (2003); (4) HST program 9055 
PI D. Charbonneau; (5) HST program 9447 PI D. Charbonneau; (6) HST program 9832 PI T. Brown. 
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Fig. 1. — Theoretical spectra of HD209458b. The upper panel shows the 1 to 6 micron window which is expected to contain many H2O 
absorption features, and the lower panel shows the spectrum from 1 to 30 microns with log flux on the y axis. Black curve: model 1, 
T eq = 1700 K, / = 2 cloud-free model; blue curve: model 2, T eq = 1450 K, / = 1 cloud-free model; red curve: model 3, T eq = 1420 K, / = 1 
cloudy model. Note how the red curve has very weak absorption bands. The dashed curves show blackbody spectra with T = 1700 K (black 
curve) and T = 1430 K (red curve); HD209458b data points (red solid error bars la, red dotted error bars 3a) are shown as follows. IRTF 
AF2.2/j m upper limit: the absorption band depth on either side of the 2.2 (im flux peak must be smaller than the vertical difference between 
the horizontal bars shown (indicated by the arrows). Because Ai<2.2(im is a relative measurement, these horizontal bars can be moved along 
the y axis. Spitzer MIPS 24/im point is shown. A hot blackbody is ruled out by the MIPS data, while model 1 (T eq = 1700 K) is ruled out 
by AF 2 . 2^im- Spitzer IRAC TrES-1 4.5/im and 8^m data are scaled for comparison with HD209458b (magenta triangles and error bars). A 
T = 1060 K blackbody (magenta dashed curve) is shown that fits the IRAC data. The renormalized Spitzer IRAC and MIPS band passes 
are shown as green dotted lines. 
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Fig. 2. — The dayside T eq as a function of Ab for different values of / (see equation 1). Top panel: T eq for the HD209458 system parameters. 
The 24 (im brightness temperature T24 = 1130 K is shown (dash-dot line); if T eq = T24 the planet must have a very high Bond albedo. 
Bottom panel: T eq for the TrES-1 system parameters. The Spitzcr IRAC measured brightness temperature of 1050 K is shown (dash-dot 
line). If T eq > 1200 K, / = 1 models are ruled out. Although the measured brightness temperatures for HD209458b and TrES-1 are similar 
(albeit measured at different wavelengths), values of T eq as a function of Ag show that the viable range of planet equilibrium temperature 
are quite different. 
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Fig. 3. — Thermal emission spectrum for HD209458b with C/O = 1.01 (and other elements in solar abundance). This chemical regime 
results in much reduced water abundance and increased CH4 opacity compared to the solar abundance spectra shown in Figure 1. The 
Ai<2.2um constraint and -F24/jm measurement are shown, as are the IRAC and MIPS bandpasses (dotted lines). 
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Fig. 4. — The photosphere and atmospheric circulation. Left panel: Pressure (as a proxy for altitude) at the planetary photosphere as a 
function of wavelength. (We define the photosphere at an optical depth of 2/3.) The f=2 cloud-free model (model 1) is shown as the solid 
curve, the f=l cloudy model (model 3) is shown as the dashed one. Right panel: the altitude dependence of the ratio of the radiative to 
advective timescales r ra( j/T at j v . A windspced U of 1000 m/s was adopted for illustration; the ratio scales linearly with U so that other values 
can be considered. See text for discussion. 



